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We study Markov chains which model genome rearrangements. 
These models are useful for studying the equilibrium distribution of 
chromosomal lengths, and are used in methods for estimating ge- 
nomic distances. The primary Markov chain studied in this paper is 
' the top-swap Markov chain. The top-swap chain is a card-shuffling 

I process with n cards divided over k decks, where the cards are ordered 

within each deck. A transition consists of choosing a random pair of 
cards, and if the cards lie in different decks, we cut each deck at the 
chosen card and exchange the tops of the two decks. We prove precise 
■ bounds on the relaxation time (inverse spectral gap) of the top-swap 

^ ' chain. In particular, we prove the relaxation time is 0{n + k). This 

resolves an open question of Durrett. 

!> 

On I 1. Introduction. Genome rearrangements play an important role in a va- 

^ ' riety of biological studies, for example, genomic distance [8, 9], phylogenetic 

. studies [11] and cancer biology [7]. Rearrangements refer to chromosomal 

fissions and fusions, reciprocal translocations between chromosomes and in- 
versions within a chromosome. Fissions, fusions and reciprocal translocations 
are examples of translocations. In this paper we study Markov chains which 
model genomic rearrangements by translocations and inversions. 

Stochastic models of chromosomal rearrangements by translocations were 
^ \ introduced by Sankoff and Feretti [10]. Such models are useful for studying 

^ ' the equilibrium distribution of chromosomal lengths. De, Durrett, Ferguson 
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and Sindi [3] studied refinements of tlie models, including the introduction 
of a fitness function on chromosomal lengths, and showed that these models 
have a reasonable fit with data for many species. These models were used by 
Durrett, Nielsen and York [6] in a Bayesian approach for estimating genomic 
distance. 

Various of the above models being simply Markov chains, Durrett [5] con- 
sidered estimating the rate of convergence to stationarity of the correspond- 
ing Markov chains. One of the open problems raised by Durrett pertains to 
the analysis of, what we call henceforth, the top-swap chain. We study the 
so-called relaxation time of the chain, which is the inverse of the spectral gap 
of the transition matrix of the chain. The relaxation time is the key quan- 
tity in the rate of convergence, and hence, it is of utmost importance in the 
efficiency of any simulations of the model. Beyond its computational signifi- 
cance, the convergence rate also has biological significance since it addresses 
the rate at which genomes reach equilibrium. 

The top-swap chain has two parameters: the number of chromosomes, 
denoted as k, and the number of genes, denoted as n. The chain can be 
viewed as a card shuffling problem with k decks. More precisely, the state 
space of the chain is a partition of n distinct cards into k decks. The cards 
have some ordering within each deck, and the decks are labeled. At each 
transition we choose two random positions, where if the size of the ith deck 
is Ui, there are + 1 positions to choose from in the ith deck. If the positions 
are in the same deck we do nothing. Otherwise, we cut both decks at the 
chosen positions. We then exchange the tops of the two decks. The figure 
below illustrates a sample transition (for k = 4,n = 12) where the chosen 
pair of positions is marked: 

4 

9 3 

4 7 2 7 5 

3 1 11 8 ^ 9 1 11 8 

5 10 6 12 2 10 6 12. 

The Markov chain allows empty decks, hence, fusions and fissions are 
modeled by this stochastic process. Translocations with nonempty decks are 
known as reciprocal translocations. As we note below, inversions can also be 
included in the above Markov chain. 

Let r(n, k) denote the inverse of the spectral gap of the top-swap chain. 
This is known as the relaxation time. We prove tight bounds, up to constant 
factors, on the spectral gap of the top-swap Markov chain. 

Theorem 1.1. There exist constants Ci > C2 > such that 
Ci{n + k)>T{n,k)>C2{n + k). 
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Remark 1.1. Durrett [5] proposes a Markov chain which is identical 
to the top-swap chain, except when the chosen pair of positions is in the 
same deck, this interval of cards is inverted. The above upper bound on 
the relaxation time immediately applies to this chain, since comparison of 
Dirichlet forms shows that extra transitions can only decrease the relaxation 
time. 



As a byproduct of our proof, we also obtain tight bounds on the relaxation 
time of a fc-deck random transposition chain; this chain is a natural extension 
of the classical (1-deck) random transposition chain to A;-decks, for k > 1. 
For a precise statement of the result, see Theorem 5.2. 

Before giving a high-level description of our proof, we need to introduce 
some notation. Let denote the state space of the top-swap chain, and 
P denote its transition matrix. We let v denote the uniform distribution 
over Q. Since P is symmetric, v is reversible. Therefore, u is the stationary 
distribution of the chain. For any / : ^ M, the Dirichlet form of the process 
is defined as 

cr,r)6f2 

and the variance is 

Var(/) = 5:K^)(/(a)-E,(/))2 

= i ^ u{aHr^){f{a)-f{rj))'. 

We then have that the relaxation time is 

T{n,k)= sup , 
/ ^U) 

where the supremum is taken over all nonconstant functions /. 

The lower bound on the relaxation time is easy to show by taking / as the 
indicator function for whether the first deck is empty. We give a high-level 
sketch of the proof outline in Section 2 before presenting a detailed proof in 
Sections 3, 4 and 5. 

2. High-level proof description. The analysis of the spectral gap of the 
top-swap chain has two major parts. The first part shows that it suffices 
to analyze the spectral gap for the 2-deck version of top-swap. (The 2- 
deck version of top-swap is simply the top-swap chain with k = 2.) We then 
analyze the 2-deck top-swap chain in the second part of our proof. Within 
both parts of the analysis, our proof relies on comparison arguments with 
other auxiliary chains. 
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2.1. Reducing to two decks. We use 2 auxiliary chains: a weighted and 
unweighted deck- averaging process. The use of a deck- averaging process is 
similar to the proof approach of Cancrini, Caputo and Martinelli [1] for the 
analysis of the so-called L-reversal chain. 

The deck-averaging process is a continuous time Markov chain where the 
state space is again all the possible partitions of the n cards into k decks. For 
each pair of decks, there is an independent Poisson clock. In the unweighted 
process, all clocks have the same rate, whereas for the weighted process, 
the clock for decks i and j has rate rij + nj , where ni denotes the number 
of cards in deck i. The proof approach of [1] immediately yields that the 
unweighted deck-averaging process has a spectral gap of 1. 

Our goal is to express the Dirichlet form and variance of the top-swap 
chain in terms of the sum of the 2-deck projections of the top-swap chain. 
It is straightforward to bound the Dirichlet form in terms of a sum over 
2-deck Dirichlet forms, that is, where the configuration is fixed outside the 
two decks. Recall that the variance is independent of the chain. Thus, we 
can bound the variance by looking at the spectral gap of the unweighted 
deck- averaging process. This leaves us with the Dirichlet form for the un- 
weighted deck- averaging process. We then use a nontrivial comparison argu- 
ment to obtain the Dirichlet form for the weighted deck- averaging process. 
Finally, it is straightforward to bound the Dirichlet form of the weighted 
deck- averaging process as a sum over 2-deck variances. The result is a bound 
on the spectral gap of the (/c-deck) top-swap chain in terms of the spectral 
gap of the 2-deck top-swap chain. 

2.2. Analysis of the 2-deck top-swap. The basic idea is to compare the 
2-deck top-swap with random transpositions. However, transpositions of two 
cards within the same deck are a problem. Roughly speaking, we can not 
efficiently "simulate" these transitions by top-swap transitions (see Remark 
3.1 below for more details). Hence, we consider a transposition chain which 
only allows transpositions of cards in different decks. Moreover, if one of the 
cards is at the top of either deck, then instead of a transposition, the chain 
does the corresponding top-swap transition. It is straightforward to compare 
this modified transposition chain with the top-swap chain. To analyze this 
modified transposition chain, we compare to a final chain, which we call 
here the balanced-swap chain. The balanced-swap chain is also defined on 
two decks and has two types of transitions: swapping and rearranging. The 
swapping transition changes the size of the decks. In particular, we choose 
a random position in a deck, and move the cards above it, to the top of 
the other deck. The rearranging transition randomizes the card order while 
maintaining the current deck size. However, we only perform the second 
transition if the decks are balanced, that is, close to the same size. (This 
is made precise with the notion of 5-balanced, for appropriately chosen < 
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6 <1.) It turns out that this later restriction is crucial for the comparison 
with the modified transposition chain. 

Finally, we analyze the gap of the balanced-swap chain by adapting an 
analytical argument given by Cancrini, Martinelli, Roberto and Toninelli in 
their work on kinetically constrained models [2]. 

In the following we first describe the solution to the 2-deck problem in 
Section 3, then proceed with the deck-averaging process in Section 4, and 
finally in Section 5 we return to the analysis of the A;-deck top-swap process. 
The matching upper bounds on the spectral gap are briefly discussed at the 
end of Section 5. 

3. Analysis of the 2-deck problem. We consider two decks Ii , I2 and 
we let ni,n2 denote the number of cards in each deck. Also, let n = ni + 
77.2 denote the total number of cards and let Qn be the set of all (n + 1)! 
arrangements of n labeled cards in the two decks. We denote by r/ a generic 
element of r2„ (a configuration) and call fi the uniform probability measure 
on r2„. 

The 2-deck top-swap Markov chain is described as follows. At each step 
we choose two positions r, s, with each position drawn at random from the 
n + 2 available positions, ni + 1 from the first deck and 712 + 1 from the 
second deck. As explained in the Introduction, if r and s belong to the same 
deck we do nothing. Otherwise, we swap the tops identified by positions r, s. 
The extra position added to each deck allows to swap an empty top. If the 
current configuration is 77, we denote by %^si] the updated configuration. 
The Dirichlet form of the Markov chain can then be written as 

1 n+2 

(3-1) ^2(/) = ^^— ^EMC?;,./-/n, 

^ ' r,s=l 

where Tr^sfiv) '■— f(Xr,s'n) for arbitrary functions / : r2„ M, and /_*[/] stands 
for expectation of a function / w.r.t. the uniform probability fi. The main 
result in this section is the following 0(n) estimate on the relaxation time 
r(n, 2) of the Markov chain described above. 

Theorem 3.1. There exists C > such that, for every n£N, 

(3.2) T{n,2)<Cn. 

Recall that the estimate (3.2) is equivalent to showing that, for every 
function /, 

(3.3) Var^(/) < Cn£2{f). 

We start by introducing a convenient notation. We add a card * to mark 
the separation between the two decks and we call (rj) the position of 

the * in rj; see Figure 1. 
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With this representation the top-swap transformations T^^s can be rewrit- 
ten by means of the transformations rj — > Tjj-r/, z, j = 1, . . . , n + 1, described 
below. 

3.1. Top-swap operators Tij. If i < x=k < j, we call TijT] the new config- 
uration obtained by the ordinary top-swap. Namely, we collect the cards on 
top of j from the second deck (including j) in a deck Dj and the cards on 
top of i from the first deck (including i) in a deck Di. We then swap their 
positions, namely, Dj goes above the position (i — 1) in the first deck and Di 
goes above the position (j — 1) in the second deck. Consider, for instance, 
the configuration rj = (3,4, *,5,2, 1,6) given in Figure 1 with = 3. If we 
pick i = I and j = 5, we obtain Ti^^t] = (2, 1, 6, *, 5, 3, 4). If we choose i = 2 
and j = 7, then T2jr] = (3, 6, *, 5, 2, 1, 4). 

If i = X* < j, then we move the deck Dj above the position (z — 1) in 
the first deck. Therefore in our example, li i = x^^ = 3 and, say, j = 6, then 
r3,6??=(3,4,l,6,*,5,2). 

Also, if i < = J, we move the deck Di above the last card in the sec- 
ond deck. Thus, in our example if, say, i = 1 and j = = 3, then Ti^^r] = 
(*,5,2,1,6,3,4). 

Finally, if i = j or 1 < i < j < x^: or 1 < x^, < i < j , we do nothing, that 
is, we define Tijrj = r/ in these cases. To finish the definition of Ti j for all 
i,j = l,...,n + l, we set Tj^i := Tij for every i < j. 

Note that with these definitions we may rewrite the Dirichlet form (3.1) 

as 

1 n+l 

(3-4) ^2(/) = ^7— ^EM(^^=./ -/)']' 

where Tjj/(r/) := f(Tijr]). Here we use the fact that the transformations T 
and T coincide when the chosen positions correspond to ordinary cards and 
are different if one of the tops to be swapped is empty. Note that the positions 
range from 1 to n + 2 in (3.1) and from 1 to n + 1 in (3.4). However, when 
one of the positions is the top of a deck in (3.1), then the second position 
must be in the other deck and there is no overcounting since in (3.4) we 
always perform the swap when one of the positions coincides with x^, . 

3 4 * 5 216 
• • • • • • • 

1 2 J 4 S 6 7 

Fig. 1. Two decks with n = 6 . The configuration rj = {3, 4,*, 5, 2, 1,6), with ui = 2,n2 = 4 
and Xt = 3. 
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3.2. Modified transpositions. A first step in the proof of Theorem 3.1 is to 
compare the process described above to a random transposition-type Markov 
chain that is defined as follows. Let Eij denote the ordinary transposition 
operator between positions i and j. Namely, {Eijrf)i = rjg for i,j ^ i and 
{Eiji])i = r]j, (Eijr])j = Next, define the modified transposition Eij, for 
i<j, by 

{Eij, iii<x^<j, 
Tij, if i = < j or i < = j, 
1, if i < j < X* or < f < j. 

Also, set Ei^i = 1 and Eij = Ej^i if z > j. Here Tjj is the top-swap operation 
defined previously. In words, Ei_j is an ordinary transposition if the two 
chosen cards (in position i and j resp.) are in different decks. If they are in 
the same deck, nothing happens. Finally, if one of the cards is the *, then a 
top-swap transformation is performed. In analogy with (3.4), we define the 
Dirichlet form of the modified random transpositions by 

1 

(3-6) ^2(/) = ^^—^EM(^M /-/)']• 

^ ^ i,j=l 

It turns out (see Lemma 3.4 below) that the comparison V2 < C£2 is easy. 
We then have to show that T>2 has a spectral gap not smaller than constant 
times n~^. This is accomplished by a further comparison with a new process 
with Dirichlet form J-sif) defined below. 

Remark 3.1. We observe that comparison with ordinary random trans- 
positions is not sufficient for our purpose. Let us call Sf^if) the Dirichlet 
form 

1 n+l 

associated to pure transpositions Note that cards are not required to 
belong to different decks. The natural normalization here, as in (3.6), would 
be 2(n + 1)^ instead of 2(n + 2)^. However, the latter has been preferred to 
make the analogy with (3.4) more apparent. Then it is simple to check that 
taking, for example, /(r/) = 1 if = (1, 2, . . . , n, *) and /(r/) = otherwise, 
the ratio (f) / £2{f) is of order n. 

3.3. The balanced-swap chain. In words, this new process is described as 
follows. Let us fix (5 G (0, 1/2). We have two independent Poisson processes a 
and r with mean 1 . At the arrival times of a we choose uniformly at random 
a position i and then update the current configuration rj by means of the 
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transformation -Ej^x* =Ti^x, [recall that denotes the position of the * 

in the configuration rj] . At the arrival times of r we look at the current value 
of x^{ifj and, if 

(3.7) 5<^<(i_5), 

then, conditioned on this value we put the system in equilibrium, that is, 
we choose the new configuration r/' uniformly at random among all the n! 
configurations with the * in the current position x^,{rj). If, instead, (3.7) 
does not hold, then we do nothing. We say that rj is (5-balanced if (3.7) 
holds. We call ps the /i-probability that this happens. Since /i is uniform, 
ps = l-25 + 0{l/n). 

We write Var^(/|x=K) for the variance of a function / : — > M w.r.t. 
the conditional probability /u(-|x=i,) obtained from ^ by conditioning on the 
position of the *. The Dirichlet form of the process described above is 
then given by 

1 n+l 

(3.8) J^sU) = ;u[x Var^(/|x,)] + — — - ^ /.[(T,,,./ - f)\ 

^ ' i=l 

where 



(3.9) xin) 



1, if ?7 is 5-balanced, 
0, otherwise. 



We will show (see Lemma 3.3) that one has the comparison < Cn2?2- The 
original problem (3.3) is then reduced to estimating from below the spectral 
gap of J-s- To this end, one can use a coupling argument, but the following 
proof seems to give a better estimate. 

Lemma 3.2. For all neN, f:^n^ we have 
(3.10) Var^(/)< 75.^5(7), 

where 'js ■= (1 - VI -Ps)~'^- 

Proof. We have to prove that the spectral gap A=k of the Dirichlet form 
Ts{-) satisfies 



(3.11) A, > 1 - v/T^. 

The generator associated to J- si') is the operator C acting on functions / as 

+ " 



I n+l 

(3.12) Cf = xil^iflx.] -f)+Af-f, Af := —— T,,xJ. 

1=1 
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Note that £ is self-adjoint and nonnegative in L?'{p) and 

(3.13) -/.[/£/] =^5(/). 

We consider the eigenvalues A and eigenfunctions fx of —C: 

(3.14) -Ch = \fx. 

Here A = corresponds to f\ = constant and any f\ in (3.14) with A 7^ 
satisfies //[/a] = 0. By definition, A* is the smallest nonzero A such that (3.14) 
holds. 

Observe that A = 1 is an eigenvalue. Indeed, it suffices to consider /(r/) = 
^p{x^{'q)) for an arbitrary 99 : {1, . . . , n + 1} — > M such that Y17=i '■Pi^) = C 
In this case, / is a function of x^^ only and, therefore, = / and the 

mean-zero property implies Af = 0, thus, (3.14) gives A = 1. 

We may assume without loss of generality that A,,, < 1, since otherwise 
there is nothing to prove in (3.11). We claim that is such that 

(3.15) /u[/aJx*]=0. 

We first need the following property. For any function (j( : ^ M, we have 

(3.16) ^i[Ag\x^]= ^i[g], 

that is, does not depend on the value of To prove (3.16), we 

write, for every j, 

I -1 1 ^ V4(7-.,g)l{.,=,}] 

n[Ag\x^ =3] = — — 2^ — — . 

Now we show that for any pair of positions 

(3.17) ^i[{Tijg)l^^^^jy\ = /x[5l|^^=i|]. 

Indeed, if 0^ := {-q G Vtn '■ x^{if) = ^}, then it is not hard to see that Tjj- : 0:^ 
0^ is a bijection: namely, if, for example, i < j, for every t] £ i}^^ and for every 
C G ^I'n, we have Ti^n+i-j+2S. G ^li^Tijrj G and 

Ti^n+i~j+2Tijr] = r], 

(3.18) 



Therefore, 



^[(?i,i5)l{x.=j}] = 7— aiTijrj) 



(n + 1)! 
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This proves (3.17). Using Ffj_[x^, = j] = and (3.17), we see that 

n+l 

1=1 

which imphes (3.16). 

We turn to the proof of (3.15). Taking the /i[-|rE^,]-expectation in (3.14), 
using (3.16) and //[/a«] = 0, we see that 

= /^[(/a, -AfxJ\x^] = ^[/A.k*]. 

Since A* < 1, this imphes (3.15). 
Using (3.15) in (3.14), we obtain 

(3.19) ^/a. = (1-A. + x)/a. =:V'/a,. 

Note that is a positive function of only. In particular, taking absolute 
values in (3.19), we have 

(3.20) fi[\AfxJ\x,]=^l^{x,)fi[\fxJ\x,]. 
Now, using \AhJ < A\fxJ and (3.16), 

fi[\AhJ\x,]<fi[A\hJ\x,]=fi[\fxJ]. 
Therefore, dividing by ^p and taking /x-expectation in (3.20), we arrive at 

(3.21) //[V'-i] > 1. 
Inequality (3.21) can be written explicitly as 

It is easily seen that this is satisfied iff — 2A* +ps < 0, or A* > 1 — ^/l — ps 
(since A* < 1 by assumption). This ends the proof of Lemma 3.2. □ 

Remark 3.2. The result in Lemma 3.2 shows that if ps = 1, that is, if 
there is no constraint and X = 1 i^^ <^ = 0), then the gap is equal to 1. It 
is, however, crucial for us that we are able to prove the spectral gap bound 
of Lemma 3.2 with the constraint X < !> since we rely on this constraint to 
perform the comparison of Lemma 3.3 below. 
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3.4. From balanced- swap moves to modified transpositions. The next 
lemma shows how to "simulate" balanced-swap moves by means of mod- 
ified transpositions. 

Lemma 3.3. For each 5 G (0,1/2), there exists Cs > such that, for 
every n and every function f on Qn, 

(3.22) Mf)<C5nV2{f). 

Proof. Let us fix an arbitrary value of x*. We want to estimate Var(/|x=K) 
by using only terms of the form fi[{Ei^i — /)^|x=f] with i < x^, < I. 

Since /i(-|x*) is nothing but the uniform probability on n! permutations 
(arrangements of n cards in two decks of given size), recalling the spectral 
gap of the ordinary random transposition chain (see, e.g., [1]) we know that 

(3.23) Var^(/|x,) < ^ E E ^K^^j/ " 

We now assume that 5n < x^, < {1 — 5)n, for 5 £ (0, 1/2) [see (3.7)]. Consider 
the case where i,j are such that i < x* and j < x* (both cards are in the 
first deck). We observe that, for any i> x^, we can write 

Ei,jf — f = Ei^iEj^iEi^if - f 

= {Eiif2 - f2) + {Ej4i - fi) + {E,,ef - /), 
where /i := Ei^ef and /2 := Ej^^Ei^f . Therefore, 
li[{E,J-ff\x,] 

<3fi[{E.^f2-f2f\x.] 

+ 3fi[iEj,efi - /i)V*] + MiEi,ef - f?\x,]. 

From the invariance of /i(-|x*) under transpositions we see that 

l^[{Ei/f2 - f2f\x*] = fi[{Ei^ef - ff\x^], 

fi[{Ej^di - fi?\x*] = l^[{Ej,ef - ff\x,]. 

Therefore, averaging over all n -|- 1 — positions i such that i > x^, we 
obtain 



l^[{E,jf-ff 



(n + l-x,)^^^ 
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Symmetrizing over the inequality above becomes 

i^m.jf - ffu.] 

Summing over all pairs i,j with i < and j < x^, we obtain 

since there are x^ — 1 positions j such that j < x^. 

Repeating the reasoning leading to (3.24) for pairs i,j with i > and 
J > X* , we have 

(3.25) 

9(n + l-x^) ^ ^ 2| 1 

^' i>xU<x, 

Inserting (3.24) and (3.25) in (3.23), we obtain 
Var(/|x,) < ^ E E /"[(^i-i/ - 

i<Xf j>x, 

9f {n + l-x,) (x, -1) ) 
+ n\ (x.-l) +(n + l-x^^rE EM(i?M/-/) 



i<x, j>Xt 

2 



where we used {n — t)/t + tj {n — t) < /t{n — t), for t G (0, n). Note that if 
x^ satisfies (3.7), then (n + 1 — x*)(x=k — 1) > (^(1 — S){n + 1)^. In conclusion, 
recalling that x forces x=k to satisfy (3.7), we see that, for any value of x^:, 

Taking yu-expectation and recalling the definition (3.8) of J^sif), we see that 
Note that Q < 2 + 9/[5(l - 5)] + 0(l/n). □ 
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3.5. From modified transpositions to top-swaps. Here we use the path- 
comparison technique of Diaconis and Saloff-Coste [4] to estimate the Dirich- 
let form P2 in terms of the Dirichlet form £2', see (3.4) and (3.6) for their 
definitions. The next lemma is based on the rather obvious fact that an 
ordinary transposition can be "simulated" by two top-swap moves. 

Lemma 3.4. For any function f , 

(3.27) P2(/) < 5£2{f). 

Proof. Let 77 be a given configuration of cards and let i,j be two given 
positions such that i < j- If i < x^:{r]) = j , then Eij = Tij, so we only need 
to take care of the case i < x^,{r]) < j. In this case Eij is an ordinary trans- 
position Eij and this can be written by way of two top-swap operations. 
Namely, suppose that r] S ri^, that is, 77 € On is such that x=k(7?) = £. Then it 
is not hard to see that 

(3.28) Eij'H = Ti+i^n+i-e+sTijT]. 

Note that n + i — i + 3<n + 2 since we are assuming i <x^, (77) = I. We adopt 
the convention that \in + i — l + 2> = n + 2, then Tij^i^n+i-i+zi changes the 
configuration ^ by moving the top above (and including) i + \ from deck 1 
to the top of deck 2. 

The identity (3.28) is easily understood with the help of a picture. For 
instance, consider the configuration r/ = (3, 4, *, 5, 2, 1, 6) of Figure 1. Ei ^rj 
(i.e., the transposition of cards labeled 3 and 1) is obtained by first applying 
Tifi to obtain Ti^rj = (1, 6, *, 5, 2, 3, 4) and then applying T2J. To see a case 
where n + i — £-|-3 = n-|-2, consider, for example, £'2,5?? (the transposition of 
cards labeled 4 and 2). This is obtained by doing first T2^5r] = (3, 2, 1, 6, *, 5, 4) 
and then applying Ts^g to this last configuration, since, by our convention, 
Ta^g is the same as T^^^. 

Using (3.28), we have, for all z < ^ < j, 

(i?„/(r?)-/(r?))2l|,.=,} 

(3.29) 

< 2(Tj+i^„+j_^+3/i - /i) l{x,=l} + '^{Tijf - f) l{x'*=£}5 

where /i(r/) := f(Tijr]) for every rj G 0^. An argument similar to that used 
in the proof of (3.17) gives 

fJ-HTi+l^n+i-e+sfl - fl)'^f{x,=£}] 

(3.30) 

= ^[(^i+l.n+i-^+s/ - /) l{x,=n+j-j+2}]- 
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From (3.30) it follows that 
i-i 

X! X! /^[(^i+l,n+i-£+3/l - /l)^l{x*=£}] 
i,j:i<j e=i+l 

(3-31) = X! X! /^[(^«+l,n+i-^+3/ - /)^l{x*=n+i-i+2}] 

i,j:i<j £=i+l 

i,j:i<j 

Inserting (3.31) in (3.29), we have shown that 

n+l ^ 

(3.32) ^'i(E^,Jf - f)%<^,<j}] < 4 E /"[(^^^Z - /)']• 

i,j=l i,j-i<j 

The identity 

1 '^+^ 

^^^^^ = h^TW ^ ^^^^^'^'^ " ^)'^i^<-<^}] 

^ ^ ij=l 

1 n+l 

+ (;^EM(T,../-/)^] 

and (3.32) then imply 

^2(/)<^— ^EMm,/-/)^] 

+ EV[(r.,x,/ - /)'] < 5£:2(/). 

n + 2 r □ 

3.6. Proof of Theorem 3.1 completed. As a consequence of Lemma 3.2, 
Lemma 3.3 and Lemma 3.4, we know that, for all n and all / : M, 

(3.33) Var^(/) < 5Csj5n£2{f). 

We can still choose the value of 5 G (0,1/2) {Cs is minimized at 5 = 1/2, 
while 75 is minimized at 5 = 0). 

Note that for S = 0.25 one has ps = 0.5 + 0{l/n), 75 = + 0(l/n) 
and Q < 50 + O(l/n). Therefore, (3.33) gives Var^(/) < Cn82{f) with C < 
875 + 0{l/n). This ends the proof of Theorem 3.1. 
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(3.34) ^^, = {f'^' 



3.7. Random transpositions with constraint. The analysis above can be 
used to obtain a spectral gap estimate for a random transposition model 
with the constraint that two cards need to be in two different decks to be 
transposed. This is the process with Dirichlet form T^Jif) as defined in (3.6), 
with the only difference that here the transformations Eij are given by 

Eij, ifi<x*<j, 

if i < j < X* or <i < j, 

with Ef- = 1 and Efj = EJ- if i > j. Namely, when i = x*, we simply trans- 
pose rather than doing the top-swap as in (3.5). 

Theorem 3.5. There exists C < oo such that, for every n € N, for every 
function f : 0,n "we have 

(3.35) Var^(/)<CnPj(/). 

Proof. We repeat the same arguments used in the proof of Theorem 
3.1. We only need to modify the definition of the Dirichlet form J'5{f) which 
we replace here with 

^ n+l 

(3.36) Fj{f) = ^[x Var^(/|x.)] + ^^-^ ^ l^m,.J - ff]. 

It is not difficult to check that all our arguments in Lemma 3.2 apply with no 
modifications to this case. Moreover, the same applies to Lemma 3.3. Note 
that here we do not need the extra comparison of Lemma 3.4. In particular, 
one has that the constant C in (3.35) satisfies C < 175 -|- 0(l/n). □ 

4. Deck-averages. We consider the following setting. There are k decks 
/i, . . . , Ifc and a total of n cards. The state space VL is the set of all 

{n + k-l)\ 

arrangements of n cards in k decks, r/ will denote the general random element 
of and u will be the uniform probability over Q. Also, we use vTjTy to denote 
the projection of the configuration ry onto the configuration of cards in deck 
h. 

4.1. Deck- averaging process. The deck-averaging process is the contin- 
uous time Markov chain obtained as follows. At each deck independently 
there is a rate-1 Poisson clock. When deck /j rings, we choose uniformly an- 
other deck Ij and the cards in the two decks /j U Ij are rearranged uniformly 
at random (among the two decks). The Dirichlet form of this process is 

k k 

(4.1) P(/) = i^^^p,,,/)2], 

i=lj=l 
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where the averaging gradient is defined by 

(4.2) = u[f\7r,r,,l^i,j]ia)-f{a) 

for i ^ j and a £ il.. We set Ai^if = 0. Here i^ij = v\-\'K^r],i ^ is the 
conditional probabihty obtained by freezing the configuration in all decks 
I I, £^ Also, note that 

(4.3) z.[(Ai,,/)2] = i.[Var,,,(/)], 

where Var^j- stands for the variance with respect to the probability Vij, with 
the convention that Varj_j(/) = 0. 

The result of Proposition 2.3 in [1] is that the spectral gap of the deck- 
average Markov chain equals 1. In particular, for every / : O — > M, we have 

(4.4) Var,(/)<P(/), 

with Vari/(/) denoting the variance of / w.r.t. ly. Strictly speaking, the proof 
of Proposition 2.3 in [1] does not apply here because the sizes of the decks 
are not fixed in our model. However, it is not hard to adapt that proof to 
show that the result (4.4) holds. This only requires small modifications, the 
point being that the spectrum of the operator /C used there satisfies the 
right bounds in our setting; see Lemma 4.13 below and the remark following 
it for more details. 



4.2. Weighted deck- averaging process. The weighted version of the deck- 
averaging chain is defined by the Dirichlet form 

k k 

(4.5) V{f) = ^ E E W(^^ + ^j) Var,,, (/)], 

where rij = nj(r/) stands for the number of cards in deck i. A possible inter- 
pretation for (4.5) is that each card is equipped with a rate-1 Poisson clock; 
when a card c rings a deck, Ij is chosen uniformly at random; if c G Ij, then 
nothing happens, while if c ^ Ij, that is, c G /j for some i ^ j, then all the 
cards in /j U Ij are rearranged uniformly at random. Apart from a factor 2, 
the Dirichlet form of this process is given by (4.5). The main result here is 
the following spectral gap estimate. 

Proposition 4.1. For every / : O ^ M, 

fik 

(4.6) Var,(/)<-P(/). 
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Proof. Since J2i=i''T'i = (4-4) and (4.3) imply 

(4.7) Var,(/) < ^ J2 E ^["^ Var,,,(/)]. 

i,j=l i=l 

We claim that, for every fixed triple i,j,i, 

^i^luiVaiijif)] < u[{ni + nj)YaTij{f)] 

(4.8) +i/[(ni + n^)Vari,K/)] 

+ i/[(nj+n£)Var J- £(/)]. 

Once we have (4.8), from the estimate in (4.7), we obtain 

Var,(/) < - E i^[ini + nj)\anj{f)], 

which is the same as (4.6). 

We turn to the proof of the claim (4.8). We can assume that are 
three distinct labels (i.e., i i, j, i/j) since the statement is obviously 
true otherwise. Let m = Ui + rij + and let ;Ua denote the conditional 
probability i/[-|7r„r/,u ^ obtained by freezing all the decks v / 
We are going to prove that 
Tn 

— Var^^ (/) < fiAiirii + Uj) Varij(/)] 

(4.9) +fiA[{ni + ne)YaTi4f)] 

+ /iA[(?^j + '^£) Var J- £(/)]. 

Here Var^^(/) is the variance of / w.r.t. /xa- This depends on the given 
configuration r] through the projections tt^t] for v ^ Note that once 

the variables iTyrj, v ^ are given, m becomes a constant, that is, m is 
//A-a.s. constant. 

Let us first show that (4.9) implies the claim (4.8). Take expectation w.r.t. 
u in (4.9) and use the fact that zv[/i] = i^[/Ua[^]] for any function /i : 17 — > M: 

ii/[mVar^^(/)] < u[{ni + nj)Yaiij{f)] 

(4.10) +i/[(ni+n^)Var,,f(/)] 

+ u[{nj +ni)Varj^e{f)]. 

Observe that 

Var^^(/)=/.A[/(/-/UA[/])] 

= fJ'Aifif - fJ-AiflT^ev])] + A^A[/(AtA[/|vr£r/] - /iA[/])] 
= /iA[Varij(/)] + Var^^ [/XAl/K^r/]] 
> ^A[Varij(/)], 
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where we use the fact that fJ-Aifl'^ev] coincides with Uij. 
We can now use the obvious bound m > to obtain 

z^[mVar^^(/)] > i/[m^A[Varij(/)]] 

= z^[?nVarjj(/)] > z^[n£ Varjj(/)]. 

With (4.10), this implies (4.8). 

It remains to prove (4.9). We introduce the projection Pi defined by 

Pif = f^A[f\T^iri]. 

Projectors Pj,Pe are defined in a similar way. Observe that with these defi- 
nitions, since ^a["K^^] we have 

(4.11) //A[Var„(/)] = ^a[/(1 - P,)/]. 
Moreover, setting 

r := (m - ?ii)(l - Pi) + - nj){l - Pj) + (m - n^)(l - Pe), 
we see that the r.h.s. of display (4.9) coincides with the quadratic form 

We call P the average: 

(4.12) p= + + 
Lemma 4.2 below shows that 

(4.13) ^tA[fPf]<h^^A[f% 

for any function / such that fJ,A[f] = 0. For such /, the result (4.13) implies 
fiAlfTf] = 2mfiA[f] - 3mfiA[fPf]+f^A[nifPif] 
+ ^J-AlrijfPjf] + lJ-A[ntfPtf] 



where we use 



L r nOn 



^lA[n^fPif]=|^A[ni{P^ff]>Q. 



The claim (4.9) follows since, by subtracting the mean /iA[/]> we may restrict 
to mean zero functions to obtain 

Var^^(/)<^A^A[/r/]. □ 

Lemma 4.2. Let P denote the operator defined in (4-12). For every 
function f such that /iA[/] =0, we have /Ua[/-P/] < 5/iA[/^]. 
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Proof. In the case of three decks, the proof of Proposition 2.2 of [1] 
[see display (27) and display (32) there] shows that 

(4.14) ^a[/(1 - P)f] > min{i, |(1 - A2)}^a[/^], 

where A2 is the largest positive eigenvalue (other than Ao = 1; see below for 
more details) of the stochastic matrix fC{a,(3) defined by 

/C(a,/3) = ^A^i?? = PlTTj-q = a]. 

That is, /C(a, /?) is the ^A-conditional probability of having the configuration 
/? in Jj, given the configuration a in Ij. 

The spectrum of /C can be computed exactly as in Lemma 2.1 of [1]. The 
only difference here is that the size of the decks is not fixed. This does not 
pose any difficulty and we can proceed in the very same way. We have to 
compute the action of /C on functions of the kind Xri • ■ ■ Xr„ , where Xr stands 
for the indicator function of the event {the card labeled r belongs to deck 
i}. For instance, a simple computation shows that for one card r one has 

^Xr = \{l-Xr)- 

Therefore, Xr — ^ is an eigenfunction with eigenvalue Ai = — ^. Similarly, for 
two distinct cards ri , r2 , one can compute 

^XriXra = i(l " Xri)(l - Xra)- 

Following the proof of Lemma 2.1 in [1], we then obtain that the spectrum 
of /C consists of the eigenvalues \i = {—\Y, for £ = 0, 1, . . . , m, where m is 
the total number of cards in the three decks. In particular, \2 = \ so that 
(4.14) implies AiA[/(l - P)f] > ^AiA[/^] or, equivalently, 

Remark 4.1. More generally, for k decks and n cards, one obtains that 
the spectrum of /C is given by 



:0,l,...,n. 



1 



k-1 



5. Top-swap process with k decks. The setting is as in the previous 
section. In analogy with the 2-deck process defined in (3.1), we consider 
here the Dirichlet form 

1 n+fc 

(5-1) W) = TTT—T^ y. ^[{%;sf - f)% 



2in + kY 
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where we sum over all positions r, s between 1 and n + k, since we added to 
each deck Ij an extra position to allow the swap of an empty top. We shall 
call Ij the set of positions defined by the deck Ij and its extra position. 
The main result is that the relaxation time of this chain is 0{n + k): 

Theorem 5.1. There exists a constant C > such that, for every n,k, 
every f 



(5.2) 



Var,(/)<C(n + A:)4(/). 



Proof. We start by decomposing into 2-deck terms: 

k 



(5.3) 



£k{f) 



1 



4(n + ky 



E i'^r,sf-f)' 



If Vij denotes the conditional probability in (4.2) (obtained by freezing all 
the decks Ii, i i,j), then 



Skif) 



(5.4) 



4:(n + kY 
1 

2{n + kf 



E iy[{ni + nj + 2f8i,,{f)l 



where £ij{f) is the Dirichlet form of the top-swap process on the two decks 
i,j only: 



2{ni + rij + 2)2 



From Theorem 3.1, we know that for all / : J7 — > M, for any pair of distinct 
decks Ii , Ij with fixed total number of cards rii + rij , 

(5.5) £i,j{f)>7Tr^, vVar,,,(/). 



C {rii + rij 



We can then use (5.4) to obtain 



(5.6) 



£k{f) > 



1 



2C{n + kf 



E y[{ni + nj + 2)\wij{f)]. 



Recalling the definition (4.5) of the weighted deck-averaging Dirichlet form 
V{f) we see that 



(5.7) 



£k{f) > 



1 



2C{n + kY 



kV{f) + 2 E ^(Var„(/)) 

i,j=l 
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Now using (4.3), (4.1) and (4.4), we get that 

(5.8) Skif) > i^c{^^W) + 2Var,(/)]. 

Proposition 4.1 then imphes 



2C{n + kf 

n+12k 
l2C{n + k)^ 



ft 

— Var,(/) + 2Var,(/) 



Var,(/) 



This concludes the proof. □ 

5.1. k-deck random transpositions with constraints. Here we consider the 
extension to k decks of the problem discussed in Theorem 3.5. Namely, we 
define 

-| n+k 

(5-9) Slif) = ^7—^ E -i(Er,sf - fn 

^ ' r,s=l 

where the transformations Er^s are interpreted as ordinary transpositions 
provided r and s belong to different extended decks If r, s G li for 

some i, then E^^s = 1- 

The computations in the proof given above can be repeated step by step. 
Using Theorem 3.5, we can then obtain the following analogue of Theorem 
5.1. 

Theorem 5.2. There exists C > such that for every n, k and every f 
(5.10) Var,(/)<C(n + fe)£:J(/). 

5.2. Matching upper bounds on spectral gap. Consider the /c-deck top- 
swap Dirichlet form given by (5.1). We now show that, up to a constant 
factor, Theorem 5.1 is tight. 

For rj £Q, let ni{rj) denote as usual the number of cards in the first deck. 
Consider the indicator (test) function 

(5-11) /(^) = (n 'l^'f^l^n 

^ ^ •'^'^ [0, ifni(?/)>0. 

Then it is easy to check that z/[/] = {k — l)/(n + A; — 1), and that 

Var,,(/) = n{k -l)/{n + k- if. 
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On the other hand, 



2{k-l)n 



2(n + fc)2 [(n + fc- 1) 

(5.12) 

n{k-l) 



(n + Zc- l)(n + /t)2' 
thus, showing that the spectral gap is at most 

{n + k-l)/{n + kf < l/{n + k). 

Note that in the first equahty above we used the fact that, for r] with ^1(77) = 
0, the contribution to the Dirichlet form comes from choosing r = 1 (viz., 
the position of the first marker) and s to be any of the n card positions (viz., 
a nonmarker position). 

Finally, note that the same test function (with an identical computation) 
shows that the spectral gap of the A;-deck random transpositions with con- 
straints, defined using the Dirichlet form (5.9), is also at most l/{n-\-k). 

Acknowledgment. We thank Fabio Martinelli for pointing out the con- 
nection with [2] and for suggesting the argument of Lemma 3.2. 
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